DOUBLE PRECISION function hitung_epsilon_mesin()
	implicit none
	DOUBLE PRECISION eps
	eps = 1
	do while ((eps+1).GT.1)
		eps = eps/2
	end do
	hitung_epsilon_mesin = eps * 2
END

DOUBLE PRECISION function contoh(x)
	implicit none
	DOUBLE PRECISION x, t
	t = 1
	contoh = exp(t)**x-5*x**2
END

DOUBLE PRECISION function contoh_aksen(x)
	implicit none
	DOUBLE PRECISION x, t
	t = 1
	contoh_aksen = exp(t)**x-10*x
END

DOUBLE PRECISION function contoh_double_aksen(x)
	implicit none
	DOUBLE PRECISION x, t
	t = 1
	contoh_double_aksen = exp(t)**x-10
END


DOUBLE PRECISION function form1A(x)
	implicit none
	DOUBLE PRECISION x
	form1A = ((x * (2.1 - 0.5 * x)**0.5) / ((1 - x) * (1.1 - 0.5 * x)**0.5)) - 3.69
END
